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I.  SUMMARY 


We  have  pursued  three  major  lines  of  investigation  this  year:  (1) 

development  of  theoretical  methods  for  predicting  the  relaxations  of  atoms  at 

IH-V  surfaces.  (2)  prediction  of  the  vibrational  spectra  of  random  III-V  and 
Hgj^Cd^Te  alloys,  and  (3)  development  of  a  theory  of  Raman  scattering  in 
correlated  clustered  alloys.  Among  other  predlcltons,  we  have  shown  that 
(1)  The  reconstruction  of  (110)  zincblende  surfaces  is  Inhibited  by  ionic 
forces,  and  so  the  29°  surface  anion  relaxation  angle  of  GaAs,  which  was 
thought  to  be  common  to  all  zincblende  (110)  surfaces,  decreases  significantly 
with  increasing  ioniclty  of  the  semiconductor.  (2)  Wur trite  (1010)  surfaces, 
to  a  reasonable  approximation,  do  not  reconstruct.  (3)  The  "clustering  mode” 
observed  in  Hg1_xCdkTe  Raman  spectra  is  due  to  Te  atoms  surrounded  by  three  Hg 
atoms  and  one  Cd  atom,  not  to  vacancy  complexes.  (4)  The  effects  of 

correlations  and  clustering  in  the  Raman  spectra  of  alloys  can  be  easily 

l 

predicted  using  a  theory  which  combines  I sing -Monte  Carlo  techniques  and  the 
recursion  method. 


II.  RECURSION  THEORY  OF  RANDOM  ALLOYS 

We  have  developed  recursion  method  calculations  for  treating  the 
electronic  and  vibrational  properties  of  alloys.  In  particular,  we  have 
predicted  the  electronic  structures  of  Pb^ _xSnxTe  alloys,  and  shown  that  the 
virtual  crystal  approximation  breaks  down  for  the  cation  s- levels,  but  not  for 
most  of  the  other  electronically  relevant  bands  (1].  To  our  knowledge,  this  is 
the  first  theory  of  this  sort  for  IV-VI  alloys,  materials  with  rocksalt 
crystal  structure  and  direct  band  gaps  at  the  L-polnt  of  the  Brlllouin  zone. 


We  have  applied  recursion  techniques  to  the  vibrational  spectra  of 
Hgi..xCdxTe  and  found  that  the  principal  features  of  the  EUutan  scattering  data 
can  be  understood,  and  that  formerly  oyster ious  features  of  the  data  can  be 
attributed  to  vibrations  of  specific  clusters  in  the  alloy.  For  example, 
Aaltharaj  et  al.  discovered  a  "clustering  node"  in  HgQ  gCdg  2Te*  vhich  they 
suspected  night  be  related  to  vacancies.  We  find  that  this  node  is  due  to  a 
vibrating  Te  atom  surrounded  by  three  Hg  atone  and  one  Cd  aton  [2].  This 
ability  to  understand  the  spectra  of  technologically  important  alloys  should 
be  useful  for  characterizing  the  materials ,  and  determining  the  importance  of 
non- random  clustering  of  the  atone  in  the  alloy. 

III.  Islng-Monte  Carlo  Theory 

We  have  developed  a  theory  of  lattice  vibrations  of  correlated  alloys. 
The  usual  assumption  one  makes  in  treating  an  alloy  such  as  Hg1.xCdxTe  is  that 
the  atoms  occupy  sites  randomly.  This  assumption  is  rarely  valid,  but  is 
invariably  made  for  lack  of  a  theory  of  non- randomness. 

Our  approach  is  to  assume  that  extended  x-ray  absorption  fine  structure 
data  (EXAFS)  that  are  nov  available  at  Notre  Dame  will  provide  information  on 
pair  correlations  in  alloys.  We  have  constructed  a  theory  which  uses  these 
data  to  construct  all  the  higher-order  correlations  using  an  Ising-Monte  Carlo 
scheme.  We  adjust  the  "exchange  integral"  and  "magnetic  field”  of  our  Ising 
model  to  yield  the  measured  pair  correlations,  and  then  predict  the 
higher-order  correlations  and  configurations  of  the  non- random  alloy.  With 
these  configurations  we  compute  the  Raman  spectra  using  the  recursion  method. 
This  combinations  of  Ising-Monte  Carlo  and  recursion  theories  yields  what  we 


Page  3 


b«llm  Is  ths  first  rsllsbls  theory  of  the  vibrational  properties  of 
non-randon  alloys. 


XV.  SURFACE  RELAXATION 

In  order  to  bettor  understand  the  aachanical  properties  of  surfaces,  we 
have  developed  a  theory  of  surface  total  energies  and  alnialred  the  surface 
energies  of  slncblende  seel  conductors  with  respect  to  rigid  rotational  bond 
relaxations.  We  find  that  the  anions  at  the  GaAs  (110)  surface  rotate  out  of 
the  surface  through  a  large  angle,  aa  Is  well-established.  We  also  find  that 
the  rotation  is  considerably  snaller  for  acre  ionic  zlncblendes  (contrary  to 
tiie  accepted  viewpoint),  and  that  tile  nore- ionic  wurtclte  (1010)  surfaces,  to 
a  good  approxlaation.  do  not  reconstruct. 
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Electronic  structure  of  Pbi^StigTo  semiconductor  alloys 
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(Received  []) 

The  electronic  structures  of  the  pseudo-binary  alloy  senlconductors 
Pb^_x8nxTe  are  analysed  using  a  tight-binding  model  with  spin-orbit 
interaction.  The  densities  of  states  end  the  band  gape  at  the  L  point  are 
eonputed  for  both  the  effective  no  die  using  the  virtual  crystal  approximation 
and  the  realistic  media  eaploylng  the  recursion  Method,  and  the  results  are 
compared.  Both  theories  exhibit  alloying  effects  such  as  band  broadening, 
energy  shifts,  and  Dismock's  band  crossing  phenomenon.  However,  significant 
deviations  from  the  virtual -crystal  approximation  are  found  for  the 
cation-derived  s-llke  deep  valence  band  states. 
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I .  Introduction 

The  narrow-gap  IV-VI  semiconductor  compounds  and  their  pseudo-binary 
alloys  have  unique  properties.  They  have  on  the  average  five  valence  electrons 
per  atom,  small  direct  band  gaps  at  the  L  point,  and  high  static  dielectric 
constants  of  order  10^.  They  often  show  a  variety  of  anomalous  thermodynamic, 
acoustic,  and  electronic  properties  [1]  [2] .  Pb^a>xSnxTe  is  an  especially 
Interesting  semiconductor  alloy  because  the  symmetry  of  valence  and  conduction 
band  edges  of  SnTe  is  reversed  compared  to  PbTe  and  other  IV-VI 
semiconductors:  The  conduction  and  the  valence  band  edges  have  Lg  and  Lg 
symmetry  respectively  in  PbTe  and  most  other  IV-VI  semiconductors,  while  the 
ordering  is  "Dlmmock  reversed"  in  SnTe  [1][2][3].  This  has  an  interesting 
consequence:  that  the  fundamental  band  gap  closes  to  zero  at  an  intermediate 
composition  x  in  Pb]_xSnxTe  (3].  This  property  of  the  fundamental  energy  band 
gap  vanishing  for  a  selected  composition  means  that  alloys  with  compositions 
near  this  composition  exhibit  small  band  gaps  that  satisfy  the  special  needs 
for  Infrared  sources  [4]  and  detectors  [5]  in  modern  technology.  Therefore  it 
is  very  important  to  understand'  the  effects  of  alloy  disorder  on  the 
electronic  structures  of  these  technologically  important  materials. 

Recently,  Spicer  et  al.  [6]  have  reported  experiments  indicating  the 
selective  breakdown  of  the  virtual  crystal  approximation  for  deep  valence 
bands  in  Hg^^Cd^e  (which  is  a  covalent  semiconductor  alloy  containing 
"light"  (Cd)  and  "heavy"  (Hg)  atoms),  and  have  identified  that  phenomenon  as 
resulting  from  the  Hg  6s  atomic  levels  being  significantly  below  the  Cd  3s 
levels.  Also,  Haas  et  al.  (7]  have  obtained  similar  disorder  effects 
theoretically  in  Hg^.xCdxTe  using  the  coherent-potential  approximation.  Davis 
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(8]  has  also  found  large  deviations  from  virtual-crystal  behavior 

2  2 

theoretically  in  Pb^.xSrxS  where  the  cations  Pb  (configuration  6s6p)  and  Sr 
(configuration  5s*)  differ  so  much  that  an  average  cation  potential  is 
meaningless . 

The  present  work  analyzes  the  effects  of  alloy  disorder  on  the  electronic 
structures  of  the  random  alloys  Fb^.xSnxTe  using  the  recursion  method  with  a 
tight-binding  model.  Pb^.xSnxTe  is  an  interesting  material  for  this  purpose 
because  its**’  constituent  semiconductor  compounds  FbTe  and  SnTe  have  very 
similar  overall  electronic  structures,  except  for  the  Dimmock  reversal  of  the 
valence  and  conduction  band  edges;  the  alloy  contains  "light"  (Sn)  and  "heavy" 
(Pb)  cations.  Moreover,  the  electronic  band  structures  of  these  materials  have 
large  spin-orbit  splittings,  and  the  fundamental  gaps  are  not  at  the  center  of 
the  Brillouin  zone,  E-5.  Indeed,  some  authors  believe  that  FbTe  and  SnTe  are 
ionic  rather  than  covalent  materials  [9] .  Therefore  the  usual  criteria  [10] 
for  the  validity  of  the  virtual  crystal  approximation  may  not  apply. 

In  Section  II,  the  tight-binding  model  for  the  parent  semiconductors  FbTe 
and  SnTe  is  discussed,  and  the  recursion  method  is  oUtllned.  In  Section  III, 
the  results  of  the  calculations  are  presented  and  discussed.  Section  IV 
summarizes  the  conclusions . 

II.  Calculations!  Procedures 


A.  Tight-binding  model 
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It  is  well  known  that  Pb^.xSnxTe  forms  a  single  phase  pseudo-binary  alloy 
over  the  entire  composition  range  x,  with  about  2%  of  lattice  constant  change 
from  PbTe  to  SnTe.  Both  compounds  crystallize  in  the  rocksalt  structure  with 
lattice  constant  6.443  A  for  PbTe  and  6.327  A  for  SnTe  (11)  at  300  °K.  The 
electronic  structures  of  PbTe  and  SnTe  (and  other  IV-VI  compounds)  have  been 
extensively  investigated  theoretically  and  experimentally  [1][2].  A  variety  of 
computational  techniques  such  as  the  relativistic  augmented  plane  wave  (APU) 
method  [12] (13] [14] ,  orthonalized  plane  wave  (OPW)  method  [15],  empirical 
pseudopotential  method  [16] [17] [18] [19] ,  and  relativistic  Green's  function  or 
Korringa-Kohn-Rostoker  method  (KKR)  [20]  have  been  used  to  calculate  the 
electronic  band  structures  of  these  materials.  More  recently,  a 
self-consistent  relativistic  APW  calculation  for  SnTe  [21]  and 
first-principles  pseudopotential  total  energy  calculation  for  the  ground  state 
properties  and  electronic  structures  of  PbTe  and  SnTe  [22]  have  been  reported. 
Although  considerable  differences  may  exist  concerning  some  details,  such  as 
the  parity  assignments  at  the  L  point  [23]  and  gap  structures  at  critical 
points  (for  example,  some  calculations  [13] [16] [21]  showed  a  "hump  structure", 
i.e.,  the  L  point  is  not  a  minimum  or  maximum  energy  p6int,  but  a  saddle  point 
in  SnTe) ,  the  general  features  of  the  various  band  structures  mentioned  above 
are  quite  similar.  Concentrating  on  this  point  and  the  fact  that  the  recursion 
method  takes  its  most  convenient  form  in  a  tight-binding  model,  we  will  use  in 
this  work  the  empirical  tight-binding  Hamiltonian  matrix  elements  of  Lent  et 
al.  [24],  which  are  obtained  by  fitting  the  eigenvalues  of  the  tight-binding 
Hamiltonian  matrix  to  the  experimental  band  gap  at  the  L  point  and  to  band 
energies  at  symmetry  points,  as  calculated  by  Herman  et  al.  [15]. 
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Since  Che  releCivisCic  corrections  Co  Che  energies  of  heavy  materials, 
particularly  chose  including  Pb,  are  significant  [25) [26],  the  Hamiltonian 
used  for  band  calculations  should  Include  these  effects.  The  relativistic 
Hamiltonian  which  produces  the  energy  band  structure  has  the  following  form 
[12] 

H  -  (p2/2m)  +  V  +  Hso  +  h272V/8m2c2‘  -  p4/8m3c2  (1) 

where  V  is  the  periodic  crystal  potential.  The  spin-orbit  term  which  may  split 
degenerate  levels  la 

Hso  “  [(AV)xp]/4m2c2  (2) 

and  the  remaining  terms  are  the  Darwin  and  mass -velocity  terms,  respectively. 

Employing  the  ideas  of  Harrison  [27],  Chadi  [28],  and  Vogl  et  al.  [29], 
the  nearest-neighbor  tight-binding  Hamiltonian  can  be  constructed 

f 

«0  "  &,o,i  E1<a  <a,i,o,R|  +  |c,i,o,R4a>  c  <c,i,a,R+3|] 

+  ^R,K',a,i,J  [|a,i,a.S>  V^j  <c,  j  ,o,R'+a| 

+  H.c. ]  +  Hso  (3) 

where  H.c.  means  Hermltlan  conjugate,  R  are  the  lattice  vectors,  i  and  j  are 
the  localized  quasl-atomlc  orbitals  for  the  cation  and  anion,  o  is  the  spin 
index  (up  or  down),  a  and  c  refer  to  the  anion  and  cation  respectively,  and  3 
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Is  the  position  of  the  cation  relative  to  the  anion  in  any  unit  cell:  d  - 
(aj/2,0,0) .  The  spin-orbit  interaction  tern  can  be  described  by  the  following 
Hamiltonian 


[|a,i,o,R>(Aa/2)Laaa<a,J,a',R| 


+  |c, i,«7,R>(Ac/2)Lc*<?c<c, J  ,ff '  .Rf] 


(4) 


As  a  basis  set,  we  used  eighteen  quasi-atomic  orbitals  localized  on  each 
atonic  site  which  are  assumed  to  be  mutually  orthonormal ized  by  the  method  of 
Lowdin  [30]:  s,  px,  py.  pz,  <^2^2,  d3z2_r2,  d^y,  dyz,  and  dzx  for  each  spin 
up  and  down  state.  The  parameters  of  this  model  are  given  in  Ref.  [24]  ,  and 
reproduce  the  experimental  band  gaps  at  the  L  point  [3]  (0.186  eV  for  PbTe  and 
0.3  eV  for  SnTe)  as  well  as  the  calculated  band  energies  of  Ref.  [IS]  at  the 
high  symmetry  points  T,  X,  and  L.  The  resulting  band  structures  are  given  in 
Ref.  [24].  In  particular,  the  Dlmmock  reversal  of  the  band  structure  from  PbTe 
to  SnTe  is  correctly  reproduced  by  the  model. 

I 

B.  Recursion  method 

In  a  random  alloy  where  the  perfect  periodicity  is  lacking,  the  wave 
vector  k  is  no  longer  a  good  quantum  number  and  the  simplifications  of  Bloch's 
theorem  do  not  apply.  Moreover,  since  the  number  of  possible  alloy 
configurations  is  infinite,  it  is  necessary  to  develop  a  theory  that  will 
correctly  predict  the  macroscopic  physical  properties  of  an  alloy  without 
having  to  perform  ensemble  averages  of  microscopic  properties  over  impossibly 
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large  numbers  of  alloy  configurations.  Various  approaches  have  been  developed 
to  overcome  these  difficulties.  Among  these,  effective  medium  theories  such  as 
the  virtual  crystal  approximation  (VGA)  [31]  and  the  coherent  potential 
approximation  (CPA)  [6]  [7]  [10]  [32]  [33]  [34]  [35]  [36]  have  been  commonly  used  for 
the  study  of  semiconductor  alloys.  In  this  work,  we  will  use  the  recursion 
method  [8] [37]  because  it  does  not  require  any  lattice  periodicity  or 
effective  medium;  and,  if  executed  to' convergence,  does  obtain  the  correct 
alloy  density  of  states  (in  contrast  with  the  effective  medium  models).  In 
this  section,  we  shall  briefly  present  the  underlying  concept  of  the  recursion 
method. 

For  the  discussion  of  electronic  structure  of  disordered  systems,  the 
local  density  of  states  (LDOS)  is  a  well-defined  and  useful  quantity;  it  can 
be  written 

n(E,?)  -  Z  |*n(?)|2  «(E-E_) 
n 

or 


-  -(1/w)  Im  G(?,?,E),  (5) 


where 

G(r,r',E)  -  Jig  <r|  (E-H+ie)  "l|r'  >  (6) 

is  the  Green's  function,  and  ^n(r)  is  an  eigenfunction  of  the  disordered 
system  with  eigenvalue  En. 
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We  choose  as  a  basis  sec  the  localized  atonic  orbitals  centered  on  the 
various  atoms  of  a  finite  cluster,  ^n(r-R) ,  where  a  runs  over  the  s,  p,  and  d 
states.  These  are  orthono realized  Lowdin  orbitals  [30]  and  R  is  the  atomic 
site.  Now  we  relabel  them  as  |1>  (i-0,  1,  2,  ....  N,  where  N+l  is  the  total 
number  of  orbitals),  and  take  as  |0>  Che  particular  orbital  (either  s,  p,  or 
d)  on  the  particular  atomic  site  (either  anion  or  cation)  we  are  interested 
in.  In  this  basis,  the  local  density 'of  states  of  the  state  |0>  becomes 

n^E)  -  -(1/w)  Im  G0  0(E)  (7) 

where  we  define 

0lfl<O  *"  <i|(E-H+iO_1|J>  (8) 

When  we  know  the  exact  eigenvalues  and  eigenstates,  we  can  calculate  the  local 
density  of  states  by  the  use  of  equation  (5) .  In  actual  calculations  for 
realistic  models  of  disordered  systems,  this  approach  requires  computations 
that  are  complicated  and  often  impractical.  Other  apprbaches  (e.g. ,  the  method 
of  moments  or  the  recursion  method)  are  needed  to  predict  the  (local)  density 
of  states  without  directly  calculating  the  eigenvalues  and  eigenvectors. 

The  recursion  method  is  one  such  scheme,  and  uses  an  analytic  expression 
to  obtain  the  density  of  states.  Since  any  Hermitian  matrix  can  be  transformed 
into  a  real  symmetric  tridiagonal  matrix  by  unitary  transformation,  the 
Hamiltonian  matrix  elements  in  the  new  transformed  basis  \v)  (*/-0,  1 . N) 


can  be  written 
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%  v* 

(j»|H|w}  -  b^+1  v-n+l 

bp  u-fi-l 

0  otherwise 


(9) 


In  other  words,  we  have  the  following  recursion  equation 


H|*0  -  bjv-l)  +  aj./}  +  b^lv+1)  (10) 


At  the  initial  Iteration,  we  choose  |0)-|0>,  and  set  bg-O.  Then  all  the 
recursion  coefficients  *v  and  b^  (*^— 0 ,  1,  ....  N)  can  be  determined  by 
recursive  use  of  equation  (10) . 

Given  the  trldiagonallzed  Hamiltonian  (equation  (9)),  the  Green's 
function  can  be  expressed  as  a  continued- fraction 


C0.0<E> 


E 


•2  - 


(ID 


where  E  has  an  infinitesimal  positive  imaginary  part.  In  principle,  Gq  q(E) 
can  be  determined  as  exactly  as  possible  by  expanding  it  to  the  N-th  level. 
However,  in  practice  the  expansion  is  cut  off  at  some  finite  level  L,  and  the 
remainder  is  neglected  here. 
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As  a  final  step,  the  continued- fraction  fora  of  Green's  function  is 
decomposed  into  partial  fractions  and  the  imaginary  part  is  taken  to  obtain 
the  local  density  of  states.  Thus  the  local  density  of  states  Is  expressed  as 
a  sum  of  Dirac  delta  functions,  and  a  relatively  accurate  Integrated  locel 
density  of  states  is  evaluated  by  integration.  Then  the  local  density  of 
states  Itself  is  determined  by  numerical  differentiation.  Details  of  the 
method  can  be  found  in  Ref.  {37] r  computer  programs  for  executing  the 
recursion  method  are  available  [ 38 ] . 

Since  we  are  Interested  in  the  macroscopic  properties  of  the  alloy,  the 

local  density  of  states  should  be  averaged  over  all  sites  in  the  alloy.  This 

can  be  done  by  the  usual  method,  l.e. ,  finding  the  local  densities  of  states 

for  each  site  and  taking  the  average  of  them.  However,  it  is  an  extremely 

time - consuming  procedure  for  a  large  cluster.  In  practice,  it  is  better  to 

choose  an  initial  state  |0)  that  is  a  sum  over  either  all  cation  or  anion 

sites  of  4a(r-R)  with  m  the  same  for  each  site  and  with  random  signs  on  the 

various  sites.  Then  the  off-diagonal  terms  <n|G|a>  that  enter  equation  (7) 

approximately  cancel  out,  due  to  the  random  signs,  giving  a  reasonable  site 

# 

average  for  the  large  cluster. 

III.  Results  and  Discussion 

We  first  calculate  the  density  of  states  for  the  perfect  crystals  PbTe 
and  SnTe  employing  the  nearest-neighbor  tight-binding  model  discussed  in  the 
previous  section.  The  results  are  shown  in  Figs.  (1)  and  (2).  The  dot-dashed 
curve  is  the  density  of  states  obtained  by  the  Lehmann-Taut  method  [39].  In 
this  method,  the  Brillouln  zone  is  decomposed  into  a  set  of  tetrahedra,  and 


\ 


1*-  * 
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th*  Lntagratlon  over  the  Brlllouln  son*  1*  evaluated  using  an  analytic 
expression.  Th*  solid  curve  is  frosi  th*  recursion  Method.  A  12x12x12  aton 
cluster  was  generated  to  simulate  the  perfect  infinite  crystal,  and  the  local 
density  of  states  for  each  orbital  was  calculated  with  th*  initial  state 
|0)-|0>  being  a  linear  superposition  of  ^B(r-E)  frou  all  th*  sites  on  a  given 
sublattica,  with  randoa  signs  (±)  on  each  sit*.  Periodic  boundary  conditions 
were  used,  and  th*  recursion  was  iterated  to  L^Sl  levels  --  that  is,  the 
continued* fraction  expansion  of  the  Green's  function  was  cut  off  at  the  b-51 
level.  Then  all  the  local  dans  1  ties  of  states  war*  sunned  to  give  the  (global) 
density  of  states. 

The  overall  agreement  between  the  two  Methods  is  very  good,  except  for 
sone  Minor  details  such  as  th*  peak  structures  and  th*  band  gap  sneering ;  the 
differences  between  the  results  of  th*  recursion  Method  and  the  Lehmann -Taut 
Mthod  are  within  the  tolerable  rang*.  The  (-function*  like  peaks  are 
associated  with  van  Hove  singularities  (40]  due  to  th*  long-range  order.  The 
nor*  or  less  smooth  peaks  in  the  upper  valence  bands  given  by  the  recursion 
Mthod  (solid  curve)  are  partly  due  to  th*  finite  site  of  the  cluster  and 

s 

partly  due  to  th*  limited  resolution  of  the  present  method  because  of  th* 
finite  cutoff  at  Lf-51.  (V*  determined  this  by  varying  th*  sis*  of  the  cluster 
and  L. )  Another  difference  is  that  while  the  Lehmann-Taut  method  clearly  shows 
the  band  gap  to  contain  no  states,  the  band  edges  are  smared  in  th*  recursion 
Mthod.  Th*  Min  reasons  for  this  are  th*  limited  resolution  of  the  Mthod  and 
th*  incomplete  cancellation  of  th*  off-diagonal  elements  of  the  Green's 
function  due  to  th*  choice  of  randomly  phased  initial  state.  The  band  edges 
can  be  sharpened  by  choosing  th*  initial  state  on  the  center  of  a  cluster  or 
by  investigating  the  spectral  density  of  states  as  will  be  discussed  later.  In 
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Fig.  (3),  Che  contribution  of  each  orbieal  to  Che  density  of  states  of  PbTe  Is 
displayed.  The  lowest  valence  band  is  predoalnantly  anion  s-llke,  and  the 
middle  valence  band  Is  cation  s-llke.  The  upper  valence  bands  have  dominant 
anion  p-like  character,  while  the  lower  conduction  bands  are  p-like  and  cation 
derived.  This  can  be  visualized  by  the  following  simple  picture.  The  Pb  atom 
has  four  valence  electrons  (6s^6p^)  with  free-atomic  orbital  energies  -12.42 
eV  and  -6.95  eV  (relative  to  vacuum)  *for  s-  and  p-  orbitals  respectively,  and 
the  Te  atom  has  six  valence  electrons  (5s^Sp^>  with  orbital  energies  -19.05  eV 
(5s)  and  -9.79  eV  (Sp)  141],  The  two  5s  electrons  of  Te,  which  have  the  lowest 
orbital  energies,  form  an  isolated  valence  band  deep  in  energy,  and  the  two  6k 
electrons  of  Pb  form  a  middle  valence  band.  The  two  6p  electrons  of  Pb  and  the 
four  5p  electrons  of  Te  Interact  with  each  other  to  form  bonding  (valence 
band)  and  antibonding  (conduction  band)  bands.  Therefore,  alloying  PbTe  and 
SnTe,  which  Is  equivalent  to  distributing  Pb  and  Sn  atoms  randomly  on  cation 
sites,  has  the  largest  effect  on  the  cation- like  middle  valence  band.  The 
characteristics  of  the  local  density  of  states  structure  In  SnTe  is  similar  to 
that  of  PbTe  (see  Fig.  (4))  the  Ss  and  5p  free-atomic  orbital  energies  of  Sn 
are  at  -12.97  eV  and  -7.21  eV,  respectively.  * 

We  generate  a  model  of  the  random  alloy  Fb^.xSnxTe  by  randomly  occupying 
cation  sites  by  either  Pb  (with  probability  1-x)  or  Sn  (with  probability  x) , 
while  all  anion  sites  are  occupied  by  Te.  The  matrix  elements  of  the  alloy 
Hamiltonian  are  derived  from  those  of  PbTe  and  SnTe  as  follows:  On  cation 
sites,  we  vise  either  PbTe  or  SnTe  matrix  elements,  depending  on  whether  the 
site  was  occupied  by  Pb  or  Sn.  On  Te  sites,  we  average  the  PbTe  and  SnTe 
matrix  elements,  weighting  the  average  in  proportion  to  the  number  of 
neighboring  Pb  and  Sn  atoms  to  the  Te.  Then  the  densities  of  states  for 
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Pb^Sn^T.  are  calculated  ualng  both  the  virtual  crystal  approximation  and  the 
recursion  method  for  a  number  of  compositions  x.  Again,  the  density  of  states 
Is  obtained  by  the  use  of  the  Lehsumn-Taut  method  In  the  virtual  crystal 
approximation,  and  a  12x12x12  atom  cluster  Is  used  In  the  recursion  method 
with  periodic  boundary  conditions.  In  order  to  avoid  sample  dependent  results, 
we  repeated  the  calculations  for  five  different  alloy  configurations  of  12^ 
atoms,  and  avaraged  the  densities  of  'states.  The  results  are  shown  In  Figs. 
(5),  (6),  (7),  and  (8).  The  solid  curves  represent  the  recursion  density  of 
states,  and  the  dot-dashed  curves  are  for  the  virtual  crystal  approximation 
(VGA)  results.  Both  the  virtual  crystal  approximation  and  the  recursion 
density  of  "states  show  the  alloying  effects  l.e. ,  energy  shifts  and  width 
changes  of  the  density  of  states  peaks.  However,  analysis  of  the  middle 
valence  band  near  -7  eV,  which  has  the  greatest  alloying  effects,  clearly 
reveals  the  differences  between  the  predictions  of  the  two  methods  --  the 
effects  of  disorder.  This  band  Is  a  doublet,  with  Its  low-  and  high-energy 
components  due  to  Pb  and  Sn  s-states,  respectively. 

Fortunately  the  results  we  find  agree  rather  well  with  what  Is  expected, 
based  on  the  Onodera-Toyosawa  theory  of  alloys  [10]  --  despite  the  fact  that 
that  theory,  to  our  knowledge,  has  not  been  applied  previously  to  alloys  with 
fundamental  band  gaps  at  the  L-point  of  the  Brillouln  zone.  The  density  of 
states  spectra  of  the  alloys  exhibit  some  features  chat  are  "persistent"  and 
others  that  are  "amalgamated”  In  the  terminology  of  Ref.  [10].  The  persistent 
features  are  associated  with  the  catlon-llke  middle  valence  bands:  the  Pb 
6s -like  and  Sn  5s -like  bands  that  retain  their  characters  In  the  alloy  because 
the  perfect-crystal  bands  do  not  overlap  In  energy.  The  remaining  bands  are 
"smalgamated"  and  tend  to  form  hybrids  of  the  PbTe  and  SnTe  bands  rather  Chan 
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exhibit  leparate  PbTe-  and  SnTa-like  bands.  This  amalgamation  occurs  because 
the  PbTe  and  SnTe  bands  overlap  in  energy,  and  hence  nix  in  the  alloy  [10]. 
Bands  that  Call  within  this  amalgamated  regime  can  generally  be  described,  in 
a  first  approxiaatlon ,  by  the  virtual  crystal  approximation. 

Although  it  is  straightforward  to  Include  a  valence  band  offset  in  the 
calculation  by  adding  a  constant  energy  to  all  of  the  diagonal  matrix  elements 
of  either  PbTe  or  SnTe  (by  construction,  the  matrix  elements  of  Ref.  [24] 
place  the  xero  of  energy  at  the  valence  band  maximum),  we  have  not  done  so 
hare  because  the  offset  is  thought  to  be  ssmll  (of  order  60  neV)  [42] ,  almost 
negligible  .on  the  scale  of  the  figures. 

It  Is  well  known  that  the  fundamental  band  gap  of  Pbj^St^Te  closes  at 
some  intermediate  composition  because  of  the  Inverted  band  structure  of  SnTe. 
He  calculated  E(Lg)-E(Lg)  of  Pb^.xSnxTe  as  a  function  of  composition  x  by 
diagonal iz ing  the  virtual -crystal  empirical  tight-binding  Hamiltonian  (closed 
circles  in  Fig.  (9)).  Also  the  corresponding  quantity  can  be  calculated  using 
the  recursion  method.  In  alloys,  the  translational  symmetry  is  broken,  thus 
the  wave  vector  E  la  not  a  good  quantum  number .  Howevefr,  we  still  can  define 
the  spectral  density  functions  analogous  to  the  perfect  crystal  by  the 
following  [8] 

A(C,E)  -  -  (1/w)  Jig  In  <b,i,o,E|C(E+i€>|b,i,oItb.  (12) 

where  |b,l,o,E>  is  a  normalised  Bloch  sum  over  all  unit  cells  of  orbital  1 
with  spin  o  on  each  atomic  site  b  (anion  or  cation).  Then  the  position  and 
broadening  of  the  peak  represents  the  energy  shift  and  damping  of  a  particular 


quaslpartlcle  a Cat*  of  energy  E  and  wav*  vector  C.  Sine*  Lg  (L£)  has  anion 
(cation)  p-llke  character,  a  Bloch  sun  of  (px+py+pr)/73  on  each  anion  (cation) 
alt*  at  th*  L  point  la  chos*n  aa  th*  Initial  stat*  |0>  for  Lg  (L£),  and  the 
apoctral  density  of  states  A(fc,E)  la  calculated.  Then  the  gap  la  defined  by 
th*  differences  In  the  peak  values  of  A(£,E),  l.e.,  E(L£)-E(Lg).  The 

theoretical  predictions  are  shown  also  In  Fig.  (9)  (closed  triangles)  In 
coaqiarlson  with  a  linear,  Interpolation  of  the  experimental  band  gaps  of  PbTe 
and  SnT*  [3]  (solid  line).  The  theoretical  uncertainty  In  E(Lg)-E(Lg)  Is 
-  ±  0.02  *V  for  0<x<l.  Th*  calculated  band  gap  Is  almost  a  linear  function  of 
composition  x  and  compares  well  with  th*  experimental  results. 

IV.  Summary 

The  electronic  structures  of  Fb^.^i^T*  alloys  Including  their  parent 

semiconductor  compounds  have  been  analysed  using  th*  tight-binding  model  with 

spin-orbit  Interaction.  Th*  densities  of  states  were  computed  for  both  th* 

effective  media  using  th*  virtual  crystal  approximation  and  th*  realistic 

media  employing  th*  recursion  method,  and  th*  results  were  compared.  As 

« 

expected,  both  theories  exhibited  alloying  effects  such  as  band  broadening  and 
energy  shifts.  However,  th*  two  methods  differed  In  their  predictions  for  the 
cation- derived  s- like  states,  which  experienced  the  greatest  alloying  effect. 
Th*  alloy  composition  dependence  of  th*  band  gap  at  th*  L  point  was  analyzed, 
and  exhibits  Dlmmock’s  band  crossing  phenomenon.  The  above  facts  show  that  th* 
recursion  method  is  a  useful  tool  for  th*  study  of  the  electronic  structure  of 
random  Ph^Si^T*.  and  In  particular  for  th*  catlon-llk*  middle  valence  band. 
However,  they  also  show  that  th*  virtual  crystal  approximation  provides  a 
remarkably  good  description  of  th*  electronically  Important  top  valence  and 
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bottom  conduction  band* .  Finally,  they  de none tr ate  that  the  Onodera-Toyozawa 
criteria  can  be  applied  to  **>].. xSnx**»  even  though  theae  alloys  have  their 
fundaaental  band  gaps  at  L:  the  cation- like  s-like  middle  valence  bands  are 
"persistent"  while  the  top  valence  band  and  lowest  conduction  band  are 

"asulgamated. " 
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FIGURE  CAPTIONS 

Fig.  (I)  The  virtual  crystal  approximation  (dot-dashed  curve)  and  the 
recursion  method  (solid  curve)  density  of  states  in  PbTe.  A  12x12x12  atom 
cluster  with  periodic  boundary  conditions  was  used  in  the  recursion  method. 
The  zero  of  energy  is  the  valence  band  maximum. 

Fig.  (2)  The  virtual  crystal  Approximation  (dot-dashed  curve)  and  the 
recursion  method  (solid  curve)  density  of  states  in  SnTe.  A  12x12x12  atom 
cluster  with  periodic  boundary  conditions  was  used  in  the  recursion  method. 

Fig.  (.3)  Local  density  of  states  for  cation  (dot-dashed  curve)  and 

anion  (solid  curve)  calculated  by  the  recursion  method  in  PbTe. 

Fig.  (4)  Local  density  of  states  for  cation  (dot-dashed  curve)  and 

anion  (solid  curve)  calculated  by  Che  recursion  method  in  SnTe. 

Fig.  (5)  The  virtual  crystal  approximation  (dot-dashed  curve)  and  the 

recursion  method  (solid  curve)  density  of  states  in  Pbp  .jSng  yTe.  A  12x12x12 

atom  cluster  with  periodic  boundary  condition  was  used  in  the  recursion 

« 

method. 

Fig.  (6)  The  virtual  crystal  approximation  (dot-dashed  curve)  and  the 
recursion  method  (solid  curve)  density  of  states  in  PbQ  4S11Q  gTe.  A  12x12x12 
atom  cluster  with  periodic  boundary  condition  was  used  in  the  recursion 


method. 
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Fig.  (7)  The  virtual  crystal  approximation  (dot-dashed  curve)  and  the 
recursion  method  (solid  curve)  density  of  states  in  Pbg  jSrig  jTe.  A  12x12x12 
atom  cluster  with  periodic  boundary  condition  was  used  in  the  recursion 
method. 

Fig.  (8)  The  virtual  crystal  approximation  (dot-dashed  curve)  and  the 
recursion  method  (solid  curve),  density  of  states  in  Pbg  gSng  ^Te.  A  12x12x12 
atom  cluster  with  periodic  boundary  condition  was  used  in  the  recursion 
method.  ~ 

Fig.  (9)  The  band  gap  E(Lg*)-E(Lg+)  of  Pb^.xSnxTe  versus  composition  x. 
The  closed  circles  (triangles)  are  obtained  using  the  virtual  crystal 
approximation  (the  recursion  method),  and  the  solid  line  represents  the 
interpolation  of  PbTe  and  SnTe  experimental  results  of  Ref.  [3], 
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"Clustering  nodes"  in  the  vibrational  spectra  of  Hg^^Cd^Te  alloys 

Zhuo-Vu  Fu  and  John  D.  Dow 
Department  of  Physics,  University  of  Notre  Dame 
Notre  Dame,  Indiana  46556  U.S.A. 

Abstract 

The  spectral  densities  of  states  for  phonons  in  Hgj^CdL^Te  alloys  are 
predicted  using  the  recursion  method  and  are  compared  with  Raman  scattering 
data.  The  "clustering  mode"  observed  by  Amlrtharaj  et  al.  is  shown  to  be  an 
"alloy  mode"  associated  with  vibrations  of  a  Te  atom  bonded  to  one  Cd  atom  and 
three  Hg  atoms.  Hg-^_xCdxTe  is  shown  to  be  a  "two-mode"  alloy  vibrationally . 
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Recent  Raman  scattering  data  [1]  for  HgQ  gCdg  2^e  substitutional 
crystalline  alloys  revealed  a  vibrational  mode  at  135  cm*^  that  is  neither 
Hg-like  nor  Cd-like.  This  "clustering  mode"  was  thought  to  be  caused  by 
either  non-random  clusters  in  the  alloy  or  by  vacancies  or  vacancy - complexes . 
Motivated  by  these  data,  we  have  executed  calculations  of  the  densities  of 
vibrational  states  of  Hgi^Cd^Te  random  substitutional  alloys,  using  the 
recursion  method  [2] .  As  a  result:  (i)  We  have  determined  that  the 
"clustering  mode"  [1]  is  attributable  to  clusters  of  one  Cd  atom  and  three  Hg 
atoms  bonded  to  the  same  Te;  there  is  no  need  to  invoke  vacancies  or 
non- random  clusters  to  describe  the  observed  mode;  and  (il)  We  have  found  that 
Hg^^Cd^Te  is  a  "two-mode"  alloy  (3]  in  that  the  optical  phonons  of  HgTe  and 
CdTe  persist  for  all  x. 

The  calculations  employ  the  recursion  method  [2],  using  a  model  with  only 
short-ranged  forces  [4] [5] [6].  The  method  and  model  are  essentially  the  same 
as  those  used  for  AiJCGa^_xAs  alloys  [6]  and  (GaAs)^_xGe2x  alloys  [7],  and  so 
we  do  not  repeat  the  formalism  here.  The  recursion  method  is  a  Green's 
function  theory,  and,  as  such,  has  elements  in  common  with  the  work  of  Talwar 
and  Vandevyver  [8].  The  force  constants  used  for  HgTe  and  CdTe  were  fit  to 
neutron- scattering  data  {9] [10],  and  are  given  in  Table  I.  They  produce  good 
phonon  dispersion  curves  and  vibrational  densities  of  states  [11] [12], 

Here  we  report  the  recursion  calculations  for  HgQ  gCdQ  2Te.  for  which 
Raman  scattering  data  are  available  [1].  Principal  peaks  of  the  spectra  of 
HgTe  and  CdTe  are  identified  in  Fig.  (1),  using  standard  notation.  In 
addition,  "alloy  modes,"  which  appear  in  the  spectra  of  neither  HgTe  nor  CdTe, 
are  indicated  on  Fig.  (1).  These  modes  are  associated,  in  general,  with 
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vibrations  of  specific  clusters  of  atoms  in  Hgj^Cc^Te.  Such  identifications 
of  peaks  in  the  spectra  with  particular  clusters  of  atoms  are  obtained  by 
choosing  various  specific  mini-clusters  [6] [7]  and  by  computing  the  local 
densities  of  states  for  one  atom  at  the  center  of  such  a  mini-cluster. 

Comparing  our  theoretical  results  with  the  Raman  scattering  data  of 
Amirtharaj  et  al.  [1]  (Fig.  (l))f  we  find  good  agreement  with  the  Cd-like 
optic  mode  at  156  cm~^  and  the  two  Hg*like  modes  at  122  and  144  cm*^ 
respectively.  The  experimental  "clustering  mode"  is  simply  an  "alloy  mode" 
[6]  due  to  Te  atoms  bonded  to  three  Hg  atoms  and  one  Cd  atom.  There  is  no 
need  to  invoke  non- randomness  in  the  alloy  or  to  postulate  the  existence  of 
vacancies  in  high  concentrations. 

tfe  find  that  the  optic  modes  of  HgTe  and  CdTe  persist  for  all  alloy 
compositions  x  in  Hgj^Cd^Te.  This  is  to  be  expected  on  the  basis  of  the 
criterion  of  Onodera  and  Toyozawa  [13]:  because  the  masses  of  Hg  and  Cd  are 
so  different,  the  spectrum  of  Hg^^Cd^Te  exhibits  "two-mode"  persistent 
behavior,  and  the  optical  phonon  bands  of  HgTe  and  CdTe  persist  for  all 
compositions  x  in  the  alloy  Hgj^Cd^Te.  Note  that  the  model  on  which  the 
recursion  calculations  are  based  contains  no  long-ranged  forces;  hence  the 
longitudinal  and  tranverse  optic  phonons  are  degenerate  in  this  model  at  £-0, 
the  Brillouin  zone  center.  Since  the  experimental  longitudinal -transverse 
optic  phonon  splitting  at  k-(5  is  small,  the  omission  of  the  long-ranged 
Coulomb  forces  responsible  for  the  splitting  does  not  alter  the  essential 
physics  or  impede  the  identification  of  the  principal  spectral  features. 
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In  summary,  recurs ion- method  calculations  of  the  densities  of  states  of 
Hg^^Cc^Te  alloys  provide  a  natural  Identification  of  the  "clustering  mode" 
with  vibrations  of  a  complex  involving  three  Hg  atoms  and  one  Cd  atom  bonded 
to  the  same  Te  atom.  A  complete  discussion  of  phonons  in  these  interesting 
alloys  will  be  published  soon  [11]. 
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£s  b  i  l<n,b,i|s>|  5(w-«s),  where  n  labels  the  central  unit  cell  of  a 

1000-atom  cluster,  b  denotes  either  a  (anion)  or  c  (cation)  in  that  cell,  i 
runs  over  x,  y,  and  z,  and  s  labels  the  normal  modes  of  the  cluster:  us  are 
the  eigenfrequencies  and  | s>  are  the  eigenvectors  of  the  cluster.  We  used  51 
levels  of  recursion,  setting  the  fifty-first  continued- fraction  coefficient  to 
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TABLE 


Table 

units  of  10 


Force -constant  parameters  (in  units  of  dyne/cm)  and  masses  (in 
^g)  for  HgTe  and  CdTe,  in  the  notation  of  Ref.  [6]. 
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FIGURE  CAPTIONS 

Fig.  (1).  The  Raman  scattering  data  of  Hgj^Cd^Te,  according  to  Ref. 

[1]  (upper  panel),  in  comparison  with  the  calculated  density  of  states  of 

H«0.8Cd0.2Te  (l°wer  panel).  The  data  were  taken  for  the  <110>  face  at  77K  and 
are  the  curves  2a,  2e,  and  2c  (from  the  top  down)  of  Fig.  2  from  Ref.  [1]. 

They  correspond  to  the  following  experimental  conditions:  X2tY2,Y2]X2, 

X2 [ Z2 *  ,Z2*  1*2*  and  XglYj.Zj]^,  where  we  have  X2[UO],  Y2[110],  Z2[001],  and 
Z2'[112].  Principal  spectral  features  are  identified  on  the  figure. 
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Abstract 

It  is  argued  that  the  surface  relaxation  angle  w  of  zlncblende  (110) 
surfaces  should  depend  on  ionicity  or  on  longitudinal  effective  charge  Z, 

0  0  o 

approximately  as  w  -  uq  -  Z^eveaL*  with  m  6  /eV. 
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The  (110)  zincblende  surface  relaxations  are  thought  to  be  the  best 
understood  of  the  sealconductor  surface  reconstructions ,  based  on  analyses  of 
numerous  low-energy  electron  diffraction  (LEED)  data 

[I]  [2] [3] [4] [5] (6] [7] (8] [9] [10]  and  medium- energy  ion-blocking  experiments 

[II]  [12].  It  Is  widely  believed  that  these  reconstructions  are  determined 
virtually  exclusively  by  covalent  forces,  and  that  the  surface  anions  rotate 
more -or -less  rigidly  out  of  the  surface  through  an  angle  u>  *■  29°  for  all 
zincblende  (110)  surfaces  [13] [14] [1] [2] .  (See  Fig.  1.)  This  relaxation  angle 
w  minimizes  the  energy  of  the  covalent  bonds.  To  be  sure,  some  (110) 
zincblende  surfaces  have  relaxations  that  are  not  purely  rigid  rotations  of 
anions  through  29°:  (i)  an  opposite  smaller  rotation  of  the  second  layer 
followed  by  some  inward  relaxation  of  the  top-bilayer  toward  the  substrate 
[10]  has  been  reported  for  some  compound  semiconductors  [9] [10] [15];  and  (11) 
small  variations  In  w  for  different  zlncblendes  have  been  purportedly 
correlated  with  variations  of  bulk  lattice  constants  rather  than  loniclty  (as 
previously  proposed  [16])  [18].  Nevertheless,  the  currently  most-accepted 
viewpoint  Is  that  the  (110)  zincblende  surface  relaxations  are  remarkably  well 
described  by  a  rigid  rotation  model  with  a  relaxation  angle  u  very  near  to  29° 
and  virtually  independent  of  ionicity. 

However,  the  viewpoint  that  w  is  constant  for  all  zlncblendes  cannot 
possibly  be  precisely  correct,  and  must  be  strictly  valid  only  in  the  limit  of 
low  ionicity.  For  the  more  ionic  zincblende  semiconductors,  the  attractive 
Coulomb  forces  between  the  relaxed,  negatively  charged  surface  anions  and  the 
positively  charged  cations  in  and  below  the  surface  tend  to  pull  the  anions 
back  toward  the  surface  and  reduce  the  relaxation  caused  by  the  covalent 
forces.  (See  Fig.  1.)  This  Coulombic  effect  should  reduce  the  relaxation  angle 


Page  3 


u  by  an  amount  proportional  to  Z2,  where  Z  is  the  (longitudinal)  effective 

2  2 

charge  (19],  and  the  effect  should  become  significant  when  Z  e  /e  a^  is  of 
order  O.S  eV,  a  typical  covalent  energy  gained  by  each  relaxing  surface  anion 
[20].  Here  e  is  the  electron's  charge,  c  is  the  dielectric  constant,  and  aL  is 
the  lattice  constant.  Recall  that  the  ionicity  f^  is  approximately 
proportional  to  Z:  Z/Zc  “  f^,  where  Zc  is  the  chemical  valency  [19]. 

Thus  it  is  disturbing  that  LEED  analyses  of  some  zlncblende 
semiconductors  with  quite  different  ioniclties  have  yielded  essentially  the 
same  values  of  u  [3] [4] [6],  especially  when  one  recognizes  that  the  electronic 
properties  of  these  surfaces  depend  dramatically  on  ionicity  [21].  Clearly  the 
surface  physics  of  zlncblende  materials  would  advance  if  either  (1)  there  were 
a  convincing  theoretical  demonstration  that  w  should  not  depend  significantly 
on  ionicity,  or  (il)  there  were  revisions  of  the  experimental  values  of  u  that 
demonstrated  a  variation  of  w  with  ionicity  or  effective  charge  Z. 

To  determine  theoretically  if  the  relaxation  angles  u  should  Indeed  be 

the  same  for  all  zlncblende  semiconductors,  we  have  executed  self-consistent 

total-energy  calculations  based  on  the  pseudo- function  local -density  theory 

[22].  These  calculations  employed  four-atom- thick  slabs  of  GaAs,  InP,  ZnTe, 

and  ZnS  [23],  and  allowed  their  (110)  surfaces  to  relax  via  rigid  rotations 

through  angles  u  until  the  total  energies  of  the  slabs  reached  minima.  The 

resulting  equilibrium  relaxation  angles  u  are  plotted  in  Fig.  2  as  functions 
2  2 

of  Z  «/(«[_  (See  Table  I.),  and,  to  a  good  approximation,  form  a  straight  line 
described  by  the  empirical  rule 

w  -  w0  -  u>1  Z2e2/«*L. 


s< 


with  Mq  -  20.62°  and  -  6.08°/eV.  Tha  fact  that  wq  la  of  the  Sana  general 
Magnitude  aa  tha  obaerved  aurfaca  relaxation  angles  indicates  that  the  rigid 
rotation  nodal,  as  used  here,  contains  tha  aaaential  physics  of  the  surface 
relaxation.  We  note  also  that  tha  values  w  we  find  for  GaAs  and  other 
zincblendes  are  conparable  with  other  theoretical  values  [9] [17] [24] [20] (25] . 
Furthermore ,  the  sign  and  magnitude  of  are  as  expected.  While  our  computed 
relaxation  angles  do  vary  systematically  with  effective  charge  or  lonlclty, 
they  show  no  such  trend  with  lattice  constant  exclusively  (See  Table  I)  and 
therefore  cast  doubt  on  data  analyses  claiming  to  have  determined  such  a 
correlations  [18]. 

The  relaxation  angles  w  extracted  from  data  analyses  are  difficult  to 
obtain  accurately  [10]  (See  Table  I.),  and  actually  agree  about  equally  well 
with  either  hypothesis;  (1)  the  widely  believed  w  -  constant  or  (2) 
k  *  «0  •  Z2e2/«aL,  with  of  order  6°/eV.  (See  Fig.  2.)  Therefore  we  think 
that  the  u  ■  constant  hypothesis  should  be  reexamined  experimentally,  possibly 
by  studying  the  (110)  surfaces  of  zlncblande  CuCi,  ZnS,  Cul,  and  ZnSe  (or  some 
other  highly  ionic  zlncblende  materials)  in  the  light  of  the  present  results. 

The  lonlclty  of  the  zlncblende  surface  not  only  reduces  the  equilibrium 
relaxation  angle  w,  but  it  also  increases  the  relaxation  energy  per  surface 
atom  and  increases  surface  phonon  frequencies,  according  to  the  theory.  (See 
Fig.  3.) 

We  hope  that  experiments  on  the  (110)  surfaces  of  different  zlncblende 
semiconductors  will  confirm  these  theoretical  ideas  quantitatively,  establish 
the  role  of  effective  charge  and  lonlclty  in  determining  equilibrium  surface 
geometry,  and  lead  to  an  even  better  understanding  of  these  prototypical 


surfaces . 
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the  bulk  that  pull  the  surface  anions  back  toward  the  bulk  through  the  Coulomb 
force  are  not  included  in  the  model.  Not  surprisingly,  both  Ref.  (24)  and  Ref. 
(25)  find  larger  values  of  w  than  we  do. 


Table  I.  Effective  bulk  (longitudinal)  charges  Z  [19),  lattice  constants 

2  2 

aL  (in  units  of  Bohr  radii  Bq),  dielectric  constants  e  (19),  Z  e  /t aL  (in  eV) , 
calculated  relaxation  angle  ‘*rheory.  predicted  change  in  relaxation  angle  Aw 
(based  on  the  empirical  rule) ,  and  relaxation  angles  w  measured  by  low  energy 
electron  diffraction  (wj  PFn)  and  ion  back- scattering  (w£on) .  respectively.  See 


Ref.  (19).  Note  that  our  results  for  theory  are  g«ner*lly  smaller  than  those 


of  Ref. 

(24) 

and  Ref 

•  (25), 

especially  for  large 

effective  charge  Z,  because 

of  superior  features 

of  our 

model  (26). 

Materials 

Z 

aL<a0> 

€  Z2e2/«*L  ^Theory 

-  Aw 

"LEED  wion 

CuCi 

1.505 

10.215 

3.6 

1.67 

10.15° 

CuBr 

1.35 

10.753 

4.4 

1.048 

6.37° 

Cul 

1.715 

11.419 

5.2 

1.347 

8.19° 

Agl 

1.38 

12.232 

4.9 

0.864 

5.25° 

ZnS 

1.795 

10.222 

5.1 

1.681  11° 

10.22° 

26°  (3) 

ZnSe 

1.625 

10.710 

5.9 

1.137 

6.91° 

29°  (4) 

ZnTe 

1.66 

11.529 

7.3 

0.891  14° 

5.42° 

28°±2°  [5] 

CdTe 

1.715 

12.240 

7.3 

0.895 

5.44° 

HgTe 

2.025 

12.208 

14.0 

0.653 

3.97° 

BN 

0.78* 

6.831 

4.5 

0.538 

3.27° 

AiP 

0.93* 

10.301 

7.6 

0.301 

1.83° 

27 . 5°±3°  [6] 

AiAs 

0.81* 

10.700 

9.0 

0.185 

1.12° 

AiSb 

1.045 

11.594 

10.2 

0.251 

1.53° 

GaP 

1.185 

10.300 

8.5 

0.436 

2.65° 

27.5°  (3) 

GaAs 

1.03 

10.683 

10.9 

0.248  20° 

1.51° 

27°-31° [ 2 ]  29°±3°  (11) 

GaSb 

0.69 

11.519 

14.4 

0.078 

0.47° 

30°±2°  (5)  28 . 5°±2 . 6°  (12) 

InP 

1.485 

11.090 

9.6 

0.563  17° 

3.42° 

28.1°  (8) 

InAs 

1.085 

11.449 

12.3 

0.227 

1.38° 

31°±3°  (7)  30°±2 .4°  (12) 

InSb 

0.72 

12.243 

15.6 

0.074 

0.45° 

28.8°  (8) 

SIC 

0.775 

8.217 

6.7 

0.297 

1.81° 

*  Estimated  from  the  ionicity  f^  by  Z  -  ZCE^ 
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FIGURE  CAPTIONS 

Fig.  1.  The  rigid  rotation  model  of  the  zincblende  (110)  surface 
geometry,  as  seen  from  the  side,  with  the  surface  anions  having  relaxed 
through  a  relaxation  angle  u  —  29°.  Note  that  the  anions  rotate  up  out  of  the 
surface,  and  that  the  bulk  semiconductor  lies  in  the  lower  part  of  the  figure. 
Note  also  that  the  Coulomb  forces  between  the  negativ  y  chargeu  anions  that 
have  rotated  out  of  the  surface  and  the  positively  charged  cations  near  and 
below  them  pull  the  surface  anions  back  toward  the  surface. 

Fig.  2.  Calculated  change  Au  -  u  •  ug  of  the  relaxation  angle  u  as  a 
2  2 

function  of  Z  e  /«aL  in  eV.  The  calculations  were  performed  for  GaAs,  InP, 

ZnTe,  and  ZnS  (closed  circles)  and  interpolated  according  to  the  empirical 

2  2 

rule  u  -  u>q  -  Z  e  /«a^  (open  circles) .  Data  for  various  semiconductors  are 
plotted  for  comparison  (See  Table  I.). 

Fig.  3.  Calculated  total  energy  curves  as  a  function  of  rotational 
relaxation  angle  u  for  GaAs,  ZnTe,  and  ZnS.  Note  that  the  effective  charge  Z 
Increases  from  GaAs  to  ZnTe  to  ZnS .  Note  that  the  angle  u  corresponding  to 
equilibrium  decreases  with  Increasing  ionicity  or  effective  charge  Z.  (See 
Table  I.)  Note  also  that  as  the  total  energy  curve  becomes  sharper  and  deeper, 
surface  phonon  frequencies  increase. 


0*0 


f 1 lefdow. nanus ]wur208.rno  10  April  1987  15:47:27 

Reconstruction  of  the  non-polar  (1010)  surfaces  of  wurtzite  AiN  and  ZnS 

H.-H.  Tsai 

Department  of  Physics,  University  of  Notre  Dame 
Notre  Dame,  Indiana  46556,  U.S.A. 

R.  V.  Kasowski 

E.  I.  du  Pont  de  Nemours  and  Company 
Central  Research  and  Development  Department,  Experimental  Station 
Wilmington,  Delaware  19898,  U.S.A. 

and 

John  D.  Dow 

Department  of  Physics,  University  of  Notre  Dame 
Notre  Dame,  Indiana  46556,  U.S.A.'* 
and 

Department  of  Chemical  Engineering  and  Materials  Science 
University  of  Minnesota,  Minneapolis,  Minnesota  55455  U.S.A. 


Abstract 


Page  2 


Using  the  pseudo* function  method  and  local  density  theory,  we  minimize 
the  total  energy  of  the  (1010)  surfaces  of  AiN  and  ZnS,  within  the  context  of 
a  Rigid  Rotation  Model.  We  predict  that  the  surface  anions  rotate  inward 
toward  the  bulk  through  very  small  angles  (in  contrast  to  the  anions  of  the 
(110)  surface  of  GaAs ,  which  rotate  outward  through  large  angles  (w“29°) . 
Thus,  to  an  adequate  approximation,  the  (1010)  surfaces  of  these  relatively 
ionic  wurtzite  semiconductors  do  not  reconstruct. 

PACS  Numbers:  68.20.+t;  73.20.-r 
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In  this  communication  we  present  calculations  of  the  surface  relaxations 
of  the  non-polar  (1010)  surfaces  of  wurtzite  AiN  and  ZnS.  We  assume  a  Rigid 
Rotation  Model  of  the  surface  similar  to  that  developed  for  zincblende  (110) 
surfaces  by  Tong  and  others  (1][2],  and  calculate  the  angle  of  rotation  w  that 
minimizes  the  surface's  total  energy. 

Studies  of  zincblende  (110)  surfaces  indicate  that  the  primary  chemistry 
determining  their  reconstruction  is  simple:  surface  anions  that  were  four-fold 
co-ordinated  in  the  bulk  rotate  out  of  the  surface  in  order  to  re -hybridize 
their  orbitals  and  adopt  a  three-fold  co-ordination  with  bond  angles  that  are 
near  90°  and  characteristic  of  three-fold  co-ordinated  molecules  [3].  To  an 
adequate  approximation,  bond- lengths  are  conserved  during  this  rotational 
relaxation  of  zincblende  (110)  surfaces  (4) . 

Guided  by  the  success  of  the  Rigid  Rotation  Model  for  the  non-polar  (110) 
surfaces  of  zincblende  semiconductors,  we  propose  an  analogous  model  for  the 
non-polar  (1010)  surfaces  of  wurtzites,  in  which  the  reconstruction  is 
described  by  a  single  angle  w,  which  gives  the  equilibrium  rotation  of  the 
anions'  surface -bonds ,  as  described  in  Fig.  (1).  We  then  determine  the 
equilibrium  relaxation  angle  w  by  minimizing  the  total  energy,  using  local 
density  theory  (5]  and  the  pseudo- function  method  [6].  To  implement  this 
procedure  without  using  excessive  computer  time,  we  take  the  solid  to  be  a 
slab  four  atoms  thick,  and,  for  simplicity,  relax  both  the  top  and  bottom 
layers  symmetrically.  Bond  lengths  are  constrained  to  remain  unaltered  during 
this  relaxation.  (Such  a  procedure  successfully  described  the  physics  of  (110) 
zincblende  relaxation  (7j.)  Details  of  our  calculational  procedure  may  be 
found  in  Ref.  [6] . 
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In  the  minimum  energy  configuration,  we  find  that  the  anions  and  cations 
of  the  (1010)  surfaces  of  wurtzite  AiN  and  ZnS  are  very  nearly  on  bulk  lattice 
sites;  the  magnitudes  of  the  equilibrium  relaxation  angles  w  are  small  and 
negative ,  -2.5°  and  -2.4°  for  AiN  and  ZnS,  respectively:  the  anions  move 
slightly  Inwards  toward  the  bulk  (in  contrast  to  the  well -understood  (110) 
zincblende  surface,  in  which  the  anions  relax  out  of  the  surface  through  a 
large  angle  of  “29°).  The  energies  gained  by  this  relaxation  are  0.14  eV  and 
0.09  eV  per  surface  atom  for  AiN  and  ZnS,  respectively  (much  smaller  that  the 
0.3  to  0.5  eV  relaxation  energies  of  zincblende  (110)  surfaces). 

This  resistance  of  wurtzite  (1010)  surfaces  to  surface  relaxation  can  be 
understood  in  terms  of  simple  physics:  Surface  reconstruction  is  dictated  by  a 
competition  between  covalent  forces  and  Coulomb  forces.  Atoms  that  are 
four-fold  coordinated  in  the  bulk  are  only  three-fold  co-ordinated  at  the 
surface,  and  covalent  forces  tend  to  rotate  these  atoms  out  of  the  surface 
through  a  relaxation  angle  u  into  configurations  such  that  they  can 
re-hybridize  and  form  three-fold  co-ordinated  bonds  with  bond-angles  near  90° 
(to  accommodate  the  p  ,  p  ,  and  p  bonding  orbitals  [3]).  The  Coulomb  forces, 
however,  tend  to  counteract  the  covalent  forces  and  pull  the  surface  atoms 
back  toward  the  bulk,  because  the  nearest  ion  to  a  surface  atom  has  an 
opposite  charge  to  the  surface  atom.  If  the  Coulomb  forces  are  sufficiently 
strong  in  comparison  with  the  covalent  forces,  the  relaxation  will  be 
inhibited. 

At  zincblende  (110)  surfaces,  relaxation  energies  are  due  primarily  to 
covalent  forces  and  are  typically  0.3  to  0.5  eV  (for  large  relaxation  angles 
w) ;  moreover  the  zincblende  relaxation  angles  are  Indeed  large:  on  the  (110) 


Page  5 


surface  of  GaAs,  the  anions  rotate  by  “29°  up  out  of  the  surface  so  that  their 
bond  angles  are  nearly  right  angles  [3] [8].  Coulomb  forces  for  a  zlncblende 

structure  tend  to  pull  the  surface  atoms  down  toward  the  bulk  and  are 

2  2 

characterized  by  energies  of  order  Z  e  /<d,  where  d  is  a  surface  bond  length 
[9],  Z  Is  the  longitudinal  effective  charge  [10],  e  is  the  electron's  charge, 
and  €  is  the  dielectric  constant  [10].  For  GaAs  this  Coulomb  energy  Is  small: 
0.25eV  [9].  Thus  for  the  zincblendfe  GaAs  (110)  surface,  the  covalent  forces 
dominate  the  Coulomb  forces  and  the  anions  rotate  out  of  the  surface  through  a 
large  relaxation  angle  w-29°  so  that  the  bond  angles  of  the  anion  are  nearly 
90°  [3J. 

For  the  (1010)  surfaces  of  wurtzite  AiN  and  ZnS,  the  opposite  is  the 
case:  Coulomb  forces  dominate  the  covalent  forces.  The  corresponding  [9] 
Coulomb  energy  for  a  wurtzite  surface  Is  Z  e  /eaL,  where  aL  Is  the  lattice 
constant,  and  is  an  order  of  magnitude  larger  than  the  Coulomb  energy  for 
zincblendes:  1.38  eV  and  2.38  eV  for  AiAs  and  ZnS,  respectively.  This  means 
that,  at  these  wurtzite  (1010)  surface,  the  Coulomb  forces,  not  the  covalent 
forces,  should  determine  the  reconstruction,  and  so  the  rotation  angle  u 
should  be  considerably  smaller  than  the  29°  characteristic  of  most  zlncblende 
(110)  surfaces:  to  a  reasonable  approximation,  these  surfaces  do  not 
reconstruct.  That  is,  we  have  w  near  zero. 

In  a  point-charge  model,  the  Coulomb  forces  on  the  anions  and  cations  in 
an  unrelaxed  (1010)  surface  are  the  same  and  so  do  not,  by  themselves,  lead  to 
a  rotation  of  the  anions  either  into  or  out  of  the  surface. 
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We  attribute  the  slight  inward  rotation  of  the  anions  to  the  strong 
second-nearest-neighbor  anion-anion  interactions  combined  with  Coulomb  forces 
that  are  unbalanced  at  the  surface:  The  anion-anion  interactions  (that  is,  the 
matrix  elements  of  the  crystal  potential  between  orbitals  centered  on 
different  anion  sites)  are  strong  and  attractive,  much  stronger  than  the 
cation-cation  interactions  --  presumably  because  of  the  greater  anion-anion 
covalent  bonding.  Furthermore,  movemeht  of  the  surface  anions  in  toward  the 
bulk  as  a  result  of  these  interactions  shortens  the  distance  from  a  surface 
anion  to  its  two  third-layer  anions  (See  Fig.  1),  further  enhancing  the 
anion-anion  interactions.  This  third  layer  effect  is  larger  than  the  effect  on 
the  inward  anion  movement  of  the  four  nearest-neighbor  anions  in  the  second 
layer,  two  whose  distances  to  the  surface  anion  are  lengthened  by  the 
rotation,  and  two  of  whose  distances  are  shortened  --  layers,  two  whose 
distances  to  the  surface  anion  are  lengthened  and  cancelling  out  the  effect  in 
lowest  order.  In  any  event,  the  inward  rotation  is  very  small. 

Thus,  if  we  assume  a  Rigid  Rotation  Model  of  the  wurtzite  (1010)  surfaces 

of  AIN  and  ZnS,  we  find  that  the  relaxation  angles  u>  are  small  and  even 

» 

negative.  Since  the  model  does  not  allow  for  small  changes  of  bond  length, 
which  could  be  as  important  as  small  rotations  of  order  w*2°,  we  cannot  rule 
out  the  possibility  that  small  bond  length  changes  might  quantitatively  alter 
the  relaxation  angles.  However  we  can  conclude  that,  to  an  adequate 
approximation,  the  (1010)  surfaces  of  wurtzite  AiN  and  ZnS  do  not  reconstruct. 
Other  wurtzite  semiconductors  with  comparable  ionic ities  and  effective  charges 
should  also  have  (1010)  surfaces  that  are  comparably  well -described  by  a 


no- reconstruction  model. 


Finally  we  note  that  Duke  et  al.  (11)  have  recently  reported  calculations 
of  the  ZnSe  and  ZnS  (1010)  and  (1120)  surfaces,  and  find  results  which  are  In 
conflict  with  ours:  they  find  large  positive  bond  rotation  angles  w  and 
considerable  surface  relaxation.  We  do  not  understand  the  physical  origin  of 
their  results,  and  note  that  the  physics  we  discuss  here  Is  not  contained  In 
the  phenomenological  model  employed  by  Duke  et  al.  That  model  replaces  the 
Coulomb  forces  with  a  phenomenological  contribution  to  the  Hamiltonian  which 
depends  on  bond-length  changes.  Such  an  approximation  may  not  be  valid  at  the 
surfaces  of  ionic  compound  semiconductors.  Ue  suspect  that  their  model,  by 
omitting  the  explicit  Coulomb  forces  between  the  surface  anions  and  the 
neighboring  bulk  atoms,  also  omits  the  force  that  inhibits  the  relaxation.  As 
a  result,  their  surface  anions  rotate  unimpeded  until  they  form  bond  angles 
near  90°  --  and  achieve  a  large  relaxation  angle  w.  If  our  suspicion  is 
correct,  the  addition  of  explicit  Coulomb  forces  to  their  model  should  lead  to 
small  values  of  «. 

It  will  be  an  interesting  challenge  to  the  experimental  surface  physics 
community  to  determine  which,  If  either,  prediction  Is  correct.  Are  the  (1010) 
surface  relaxations  of  wurtzlte  AiN  and  ZnS  large  or  small? 

Acknowledgments  --  We  are  grateful  to  the  U.  S.  Air  Force  Office  of 
Scientific  Research  (Contract  No.  AF-AF0SR-85-0331)  and  E.  I.  du  Pont  de 
Nemours  and  Company  for  their  support. 


REFERENCES 


*  Permanent  address. 

(1]  S.  Y.  Tong,  A.  R.  Lublnsky,  B.  J.  Mrstik,  and  M.  A.  Van  Hove,  Phys .  Rev.  B 

17.  3303  (1978). 

(2]  D.  J.  Chadi,  Phys.  Rev.  B  18,  18000  (1978);  19,  2074  (1979). 

(3]  J.  J.  Barton,  W.  A.  Goddard  III,  and  T.  C.  McGill,  J.  Vac.  Sci.  Technol . 

16,  1173  (1979). 

(4]  P.  G.  Cowell,  M.  Prutton,  and  S.  P.  Tear,  Surf.  Sci.  177.  L915,  (1986);  A. 

Kahn,  Surf.  Sci.  168 .  1  (1986);  C.  B.  Duke,  C.  Mailhiot,  A.  Pa ton,  D.  J. 
Chadi,  and  A.  Kahn,  J.  Vac.  Sci.  Technol.  B3,  1087  (1985). 

(5]  L.  Hedin  and  B.  I.  Lundqvist,  J.  Phys.  C  4,  2064  (1971). 

(6]  R.  V.  Kasowski,  M.-H.  Tsai,  T.  N.  Rhodin,  and  D.  D.  Chambliss,  Phys.  Rev. 

B34,  2656  (1986). 

[7J  R.  V.  Kasowski,  M.-H.  Tsai,  and  J.  D.  Dow,  to  be  published. 

[8]  M.  W.  Puga,  G.  Xu,  and  S.  Y.  Tong,  Surf.  Sci.  164,  L789  (1985). 

[9]  In  terms  of  the  zincblende  latticeice  constant  a^,  we  have  d  - 

[10]  G.  Lucovsky,  R.  M.  Martin,  and  E.  Burstein,  Phys.  Rev  B4,  1367  (1971). 

[11)  C.  B.  Duke,  Y.  R.  Wang,  and  C.  Mailhiot,  Bull.  Amer.  Phys.  Soc.  32,  709 
(1987).  See  also,  C.  Mailhiot,  C.  B.  Duke,  and  D.  J.  Chadi,  Surf.  Sci. 
149,  366  (1985). 


Figure  Caption 


Fig.  (1).  Illustration  of  the  reconstruction  of  the  (1010)  wurtzite 
surface  for  a  large  positive  relaxation  angle  w.  The  anions  and  cations 
are  denoted  by  -  and  +  signs,  respectively. 


Side  view  of  the  wurtzite  C1010)  surface 


